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Abstract 

We study real-time correlation functions in scalar quantum field the- 
ories at temperature T = 1//3. We show that the behaviour of soft, 
long wave length modes is determined by classical statistical field the- 
ory. The loss of quantum coherence is due to interactions with the 
soft modes of the thermal bath. The soft modes are separated from 
the hard modes by an infrared cutoff A <C l/{fi(3). Integrating out the 
hard modes yields an effective theory for the soft modes. The infrared 
cutoff A controUs corrections to the classical limit which are O {hp A). 
As an application, the plasmon damping rate is calculated. 



1 Introduction 



The properties of quantum field theories at high temperature are important 
in order to understand the behaviour of matter under extreme conditions, 
as they occur in astrophysics and cosmology. Static quantities at high tem- 
perature characterize the thermodynamics of the system, in particular the 
equation of state. Time-dependent quantities are relevant for the study of 
non-equilibrium processes, such as the generation of a cosmological baryon 
asymmetry. 

Consider the partition function for a scalar field theory at temperature 
T = 

Z = I D^e-^T'^-'^'-^W , (1) 

where the functional integration is over periodic fields, (^{xjhp) = (/)(x, 0). 
Formally, the limit T — oo is closely related to the limit h —>■ 0. This suggests 
that the high temperature behaviour of bosonic systems may be described by 
an effective classical field theory. For static quantities this reasoning leads 
to the dimensional reduction program which has been successfully applied 
to the electroweak phase transition in recent years [|l|]. It has also been sug- 
gested that the classical theory can be used to compute real-time correlation 
functions at high temperature 0. This is of particular interest for sphaleron 
processes in the electroweak theory where, however, results for the transition 
rate are still controversial 0, ^ . 

Scalar field theories are much simpler than gauge theories and therefore 
a useful starting point to study real-time processes at high temperature [^]. 
This concerns the problem of ultraviolet divergencies p as well as the ques- 
tion of quantum corrections to the classical limit 0. In classical 0^-theory 
one studies the classical correlation functions 

(0(Xi, ti) . . . 0(X2, t2))e, = ^ / 'D7tV<P g-^^^^'-^) 0(xi, ti) . . . 0(x2, ^2) , (2) 

where 

Hin, <P) = j (^vr^ + i(V0)^ + \^^'<p' + ^A^^) (3) 

is the classical Hamiltonian, and 0(x, t) is the time-dependent field which is 
obtained by integrating the classical equations of motion from initial condi- 
tions 0(x, tj), 7r(x, ti) over which then a thermal average is performed. From 
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the quantum thermal two-point function one obtains the plasmon damping 
rate . As shown by Aarts and Smit , the damping rate is determined by 
the classical theory to leading order in the coupling A. This result is most 
easily obtained from the retarded Green function, and it turns out that the 
whole imaginary part of the self-energy is given by the classical theory in the 



high-temperature limit [T^. The classical limit for scalar field theories has 



been further discussed in 11, 12 



The quantity h(3 carries dimension. Hence, the question arises which 
dimensionless quantity controls the classical limit. As is well known, the 
spectral energy density of a relativistic Bose gas, i.e., Planck's formula, 

de = — ^ „ ^ duj , (4) 



de = ■^Tu'^duj (5) 



turns into the classical Rayleigh- Jeans formula 
for small frequencies, 

h/3iu<^l . (6) 

Hence, the classical limit is obtained for low frequency, long wave length 
modes. For high frequency modes there is no classical limit. 

We conclude that a consistent definition of the classical limit requires the 
introduction of an infrared cutoff A <^ l/{hf3) which separates 'soft' modes 
with uj < A from 'hard' modes with u > A. The effective theory for the soft 
modes, where the hard modes have been integrated out, can be approximated 
by a classical theory. For massive theories one has m < hu < HA. Since at 
finite temperature a plasmon mass ~ hXT^ is generated, the condition 
@ implies -C 1. Hence, a classical limit can only be obtained in the case 
of weak coupling. 

In the following sections we shall discuss the connection between the 
finite-temperature quantum theory and the classical theory in detail. In 
sect. 2 we construct an effective theory for the low frequency modes. This 
is similar in spirit to the construction of an effective average action [in 



and in particular to the Wilson renormalization group approach for thermal 



modes ||15[. Sects. 3 and 4 deal with the infrared approximation which yields 
the classical limit up to corrections O {h(3A) and quantum corrections. As an 
application, we discuss the evaluation of the plasmon damping rate in sect. 5, 
and sect. 6 summarizes the main results. In the following we set ^ = 1. 
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2 Low energy effective theory 



The basic quantities of finite-temperature field theory are the thermal aver- 
ages of time-ordered products of field operators, 

(t<I(xi) . . . <l(x„)) = I Tr (e-^^ T<l(xi) . . . I>(x„)) , (7) 
which can be obtained by differentiation from the functional 

Z[j] = Tr(e-^^T,e*/c'^"^*) . (8) 

Here Tc denotes ordering along an appropriately chosen path C in the com- 
plex time plane [0, and J^dxj^ = J^dt J (i^xj(x, t)<l>(x, t). 

The dynamics of modes with energies u below the temperature T is 
strongly affected by the thermal bath. In order to derive a low energy ef- 
fective action for these modes one has to integrate out hard thermal modes. 
In the following we shall perform the part of the trace in Eq. (0) which 
involves large momentum modes. This procedure is similar to the Wilson 



renormalization group approach for thermal modes [|l^ . 

A convenient basis for the evaluation of the trace is provided by the 
coherent states of Bargmann and Fock (cf. |T^). They are constructed in 
terms of creation operators related to the field operator at some fixed time 

til 

\7]) = N elk^^"'^ \0) , (9) 
where is a normalization factor and 



The coherent states are eigenstates of annihilation operators, 

a^lv) =Vk\v) ■ (11) 
In the following we shall also use coherent states restricted to soft modes, 

M = ArAe^w<A'"''^^|0) , (12) 

As discussed in the introduction, hard modes are defined by cuk > A, where 
the separaton scale A has to be chosen in the range m < A -C 
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We can now split the trace which gives the functional Z[j] into soft and 
hard modes, 



Z[j] 



Vr]Vr]*ZA[j] 



(13) 



where 



Za[j] 



VrjVrj* 1 7] 



e 2' 



oj>A 



TiA (e 



(14) 



and VrjVrj* oc nk'^'7k'^^k- The dependence of Z\[i] on the soft modes 77k, 
cuk < A, will be discussed in detail below. Note, that we have written the 
trace in a symmetric form in order to maintain the reality of expectation 
values of hermitian operators. 



ti+i/2 (3 



c. 



ti- i/2 (3 



Figure 1: Keldysh-type contour for real-time correlation functions 

The integral over the hard thermal modes yielding the functional Z/<^[j] 
can be carried out in perturbation theory. Using the usual methods of real- 
time perturbation theory one easily derives the identity 

e-^^^T,e*/c'^^^*e-^^^ = e'i'c'^'^^^(H)e-^^^« T,e*i'c'^^^*« e"^'^^" , (15) 

where the subscript refers to free fields which are related to interacting 
fields by a unitary transformation. C is the Keldysh-type contour shown in 
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Fig. g, 



From Eqs. (|I4D and (^) one reads off, 
where 



Zf[3]= I Vr^Vr^* ir, 



7] 



(16) 



(17) 



(18) 



a;>A 



In order to evaluate the integral over hard thermal modes we first express 
the time-ordered product in terms of a normal-ordered product, 



rj J^dxj^o _ ^-^ J^dxidx2j{xi)Gc(xi,X2)j{x2) . J^dxj^o 



(19) 



where Gc is the usual causal propagator |18 



(e( 



n - T2 e 



-«i^k(ii-t2) 



+0(r2-ri)e*^''(*i~*2)) g^'^^'^i-^^) . (20) 

Here r, which parametrizes the contour C, is a real variable whereas t(r) is 
in general complex. 

For matrix elements of normal ordered products one easily proves the 
identity 



^k 



This yields 



V 



-l3Ho 



(21) 
(22) 



where the field is given by 



(x,t) 



V 



0o(x,t) 



/ (0k cosci;k(t — tj) H — ^ sinci;k(i — ij) ) e*''^ 

Jk V CUk / 



(23) 
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with 



!)k = e 2 



(24) 



Here \ri) = exp (— i/5i/o)|^)- A straightforward calculation yields for the 
matrix element of the free Hamiltonian, 



V 



Tj) = exp — - 



4 Jk n[ijjy^) \ 

where n{uj) is the Bose-Einstein distribution function 

1 



Fk 



n{uj) 



(25) 



(26) 



The free functional can now be obtained from Eqs. (p!9|), ( p2D and 

(25) by performing a Gaussian integration. The calculation gives the result 

_ g-/3^fo,A g-^ /^(ia::ida:2j(xi)GA(a;i,X2)i(a;2) dxj</>A ^^27) 

where -ffo.A and 0a are the free hamiltonian and the free field, respectively, 
restricted to soft modes, 



-/3Ho,A 



0A(a;) 



^A 



(28) 



(^A 


0o(x,t) 


^a) 







/ Q{uj - A) f0kcoscjk(t - U) + — sincjk(t - U)] e^*"^ ,(29) 

Jk V CUk / 



with It^a) = exp (— l/Jffo)!^)- Integration over the hard modes leads to a 
corresponding thermal contribution to the propagator. 



(30) 



= 2 / nA(^k)cos(a;k(ti -t2))e''^^''^-''^) , (31) 

Jk 

where nA(co') = Qiu — A)n(co'). From Eqs. ([TtD and ( pTl ) one obtains 



-/^-H'o.A g* X; '^^■J"?^A 



^Ab1 , 



(32) 



6 



where 

Note, that Z\ and Z\ are functionals of j and 0a, the soft modes, which play 
a role similar to a background field. 

Using Eqs. (|14D and (|32D one can express the thermal Green functions in 
a compact form. Changing variables from rj, rj* to (p, n the result reads 



j=0 



where 



= ^ J V<pVrre~^''^{T^x,)...^Xn))^ (34) 

u}<A 

= Ho,A + Hi^A = ~ In Za[0] , (35) 



and 



(36) 



j=0 

From these expressions one can easily recover the usual real-time pertur- 
bation theory. The generating functional for the thermal Green functions is 
obtained by integrating Z\[j] over the soft modes. A Gaussian integration 
yields, 

J VrjVr]* e''^^°-^ e^-f'^'^^-''^^ ^ ^-^ J^<ixidx2j(xi)GT,K{xi-x2)j{x2) j-g-j?^ 



where 



G't,a(x) = 2 / nA(cc^k) cos(a;kt)e'''" , (38) 



where n\{uj) = n{uj) — n\{uj) = 0(A — uj)n{uj). Combining Eqs. ( p2D and 
one obtains the familiar result 

^j-jj ^ e*/c'^^^^(T^) Ld.xidx2j{xi)GixuX2)j{x2) ^ j^gg^ 



where G = G\ + Gt,a is the usual thermal propagator [|r^ . The perturbation 
series for Green functions with infrared cutoff A is obtained from the pertur- 
bation series without cutoff by the following procedure: one substitutes the 
propagator G{x — y) by Ga{x — y), and one adds contributions which are 
obtained by replacing a propagator G{x — y) by the product of background 
fields (l)A{x)(f)A{y). 
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3 Infrared approximation 



So far we have only reorganized the perturbative expansion, without making 
any approximation. In the resulting low energy effective theory, however, we 
have a small parameter, j3K <^ 1, where A is the cutoff which distinguishes 
between soft and hard modes. As discussed in the introduction, in an expan- 
sion with respect to /3A the leading term should correspond to the classical 
limit. We shall refer to the leading term as the infrared (IR) approximation 
since, as we shall see, further approximations are necessary in order to obtain 
the classical limit. 

Consider the thermal Green functions with IR cutoff A. From Eq. (^) 
one obtains 



Za[j] = e-^^°-'^ e'L"''^'^' 0a(x) + --f- Z^j] , (40) 



i 6j{x) \ i Sj{x) 

and therefore 

T<l(xi)...<l>(x„)), = ^UA(a;i) + - 



■ (41) 

j=0 



The right-hand side of the equation is the sum of products oi n — m fields 
4>A {0 < m < n) and m-point functions 



~f .lis 5 ~ 



(42) 

j=0 



The f-functions can be computed as usual in perturbation theory. Ac- 
cording to Eq. (p2D, (j)A{x) occurs as a background field which leads to ad- 
ditional vertices (cf. Fig. |^).The f-functions can be decomposed into discon- 
nected and connected parts (cf., e.g., p9[]), 

r{yu...,ym) = r{yi) . . .f{y^) 

+T{yi, y2)c riys) ■ ■ ■ r{ym) + permutations 
+ . . . 

+f{yi,...,ym)c ■ (43) 
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Figure 2: Vertices in the presence of the background field 0a 

For the tadpole one obviously has (cf. (^)), 

fiy) = (^y))^-My)- (44) 
From Eqs. (p])- (|4^) one obtains 

(tI>(xi) . . . <l(xO)^ = (Hxi))^ . . . (Hxn))^ + PC , (45) 

where PC denotes partially connected contributions which contain connected 
m-point functions. 

^ 1 \ / ^3 



Figure 3: Contribution to )r with cuts indicated 

Consider a particular contribution to PC with a connected m-point func- 
tion f{yi, . . . , ym)c (cf- Fig. Substituting for m — 1 lines 

G^{y-z)^<PM<pA{z) (46) 
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yields a contribution to . . (cf. Fig. ^). Hence, for each 

contribution to PC there exists a corresponding term contained in the dis- 



connected product . . . (^^{xr,, ^ ^ 

In order to obtain the complete thermal Green functions we still have 
to integrate over the soft modes. As discussed in the previous section, the 
integration /^^^ D^Ptt replaces products of background fields by the infrared 
propagator Gt.a (cf. Eq. 



My)M^)-^ I GTA^,yo-zo) (^''^^-'K (47) 

Ll1<A 

For soft modes, with 

cu < A < T , (48) 

the propagators G\ and Gt,a, which occur in the disconnected and connected 
contributions, respectively, satisfy the inequality, 

GA(k, t) — < pA ^ PA Gt.a{K t) . (49) 

This implies that the dominant contribution to thermal Green functions is 
determined by the disconnected term <^<l>(xi)^^ . . . ^<l>(x„)^^ for soft external 
momenta, i.e., |xj — Xj| > 1/A. From Eqs. ( pi]) and (|45| ) one obtains in this 



case, 



T<l(xi)...<l(x„)) = | J P0D7re-'^^^(<l(xi))^...(<l(x„))^ + O(/3A) . 

(50) 




Figure 4: The tadpole diagrams corresponding to the cut components 
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This equation is the main result of this paper. It is similar to the ansatz of 
classical statistical field theory, as described in the introduction. There are, 
however, several differences with respect to the classical theory. 

Consider first the hamiltonian i^A, as defined in Eqs. ( p5D and (|36D. The 
interaction part is, to leading order in the loop expansion, 

Hi,K = - j d?xCi{<P{^,U)) + 0{(3K) . (51) 

Higher order corrections can be calculated in perturbation theory and yield 
further deviations from the classical Hamiltonian. Further, as we shall see 
in the next section, also the initial conditions for the fields <^<l>(x, t)^^ and 

^n(x, t)^^ = <^(9(<l>(x, t)^^ are different from the integration variables in the 
functional integral (|50|). 

Hence, except for the momentum cutoff, the integration over (p and vr 
corresponds to an integral over all possible initial conditions. Eq. (|50| ) differs 
from the naive classical limit by the presence of the momentum cutoff A, by 
quantum corrections to the classical hamiltonian and by the deviation of the 
time evolution of ^l>(x)^^ from the classical evolution. 

It is very interesting that in the high-temperature limit the coherence of 
the quantum fields is lost to leading order in /?A. In Eq. (|^) the expectation 
value of the product of field operators has been replaced by the corresponding 
product of expectation values, and the time ordering has therefore lost its 
meaning. This is an effect of the thermal bath where the propagator of the 
soft modes is dominated by the thermal part Gt,a- 

4 Time evolution of expectation values 

In order to understand the connection between the infrared approximation 
and the classical approximation one has to study the time evolution of field 
expectation values in the low energy effective theory. The time dependent 
field 

{Hx))^ = ^TTA{e-^^Hx)) (52) 
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can be evaluated in perturbation theory. From Eqs. (|3^ ) and (^) one reads 
off, 

/ 1 5 \ 

^ix)). = Ua(s) + - 



A 



i 5j{x) J 



■ (53) 

j=0 



To be specific, consider a quartic self-interaction 

A 



CM = ^<f>\ (54) 



The corresponding Feynman diagrams contributing to (^^{x)j^ are shown in 
Fig. ^ up to two loops. 

The kinetic operator for a mode with momentum k is 

Ko{k) = ^^+ul. (55) 



From the definitions ( pQ]) and (|30D one reads off that the background field 0a 
and the propagator G\ satisfy the equations {ti < t < tf), 

iro(k)0A(k,t) = 0, 

Koik)GAik,t) = -t5it) . (56) 

Furthermore, one easily obtains for the field ^l>(x)^ , up to terms O {X'^) 
(cf. Fig. D, 



Ko ($)^ = - A - (D, + + D,) (1.)^ - ($)^ A ($)^ . (57) 

Here the Dn are non-local kernels which give the contributions of the corre- 
sponding IPI diagrams in Fig. ^ Eq. (^) is similar to the classical equation 
of motion. The only difference is that the coupling constants and A have 
been replaced by non-local kernels. Note also, that higher orders in the loop 
expansion will generate higher powers of ($) in Eq. (| 



A 

Also important are the initial conditions which solutions of Eq. ( |57D have 
to satisfy. From Eq. (|53|) it is clear that 

I'(x,O))^^0a(x,O) . (58) 
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Instead, one finds for the Fourier coefficients of both fields, 



l>(k,0))^ =0k + O(A/5A) 



dMk,t) 



7rk + 0(A/5A) . (59) 



Hence, the integration variables in the infrared approximation (^0|) agree with 



the initial conditions of the field (^^{x)^ ^ only up to corrections O {X/3A). 

Consider now the quantum corrections to the classical field equation. At 
one loop we have to consider the diagrams Figs, and The first term is 
local and leads to an effective thermal mass, 

miff = < + ^l GA(k, 0) = < + ^T'(l-^^ /3a) = m^iry . (60) 

Here rriB and are the bare mass and the renormalized mass of the zero 
temperature 4D theory, respectively. The thermal mass m\{T) is well known 
from finite-temperature perturbation theory. 
Diagram ^ yields the contribution 

Hx))^ D,{^)l{x) = -i^ («l(x))^ ldx'G^ix,xr , (61) 

where the contour integral can be carried out by means of the two-component 
formalism (cf. appendix). The divergent contribution is again renormalized 
as in the zero temperature 4D theory. The contributions §6 and |^ can be 
combined to an effective non-local coupling. A straightforward calculation 
yields the result (cf. appendix), 

^effi^ - = A(r)(5(x - X') + Xnonloc{x - x'), (62) 

where A(T) is the temperature dependent coupling familiar from dimensional 
reduction 

A(r) = XM (i - 1^ f ) + ^ (f^^) ' (63) 

and 
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^nonloc 



(k,t-t') 



3A2 
4 



(27r)3 CJqWk.q 

1 + nA(cUq) + riAiuk-d) 



+ 



cos(u;q + u;k-q)(^ - t') + 

COs(cUq-CJk-q)(t-t') e(t-t'^ 



d_ 
dt' 
(64) 



To obtain this result we have performed a partial integration in t' and dis- 
carded the contribution at t' = —oo (cf. [^). From the upper limit of 
integration t' = t one obtains a contribution which is local in space up to 
corrections O (jSA) The important infrared contribution to the non- 

local part of the coupling was previously found by Boyanovski et al. in a 
non-equilibrium context |20|. 



5 Retarded Green function 

An interesting illustration of the classical limit is the calculation of the plas- 
mon damping rate 0, which is most easily evaluated by considering the 
retarded Green function |T0| . 

The retarded Green function of the low energy effective theory is given 

by 

^D^{x^,X2) = i TrA (e"^^ [^x,),^X2)]) Qih - t^) . (65) 

It can be obtained from the time-ordered Green function by means of the 
following relations, 

tD^ixu X2) = {D>{xu X2) - D<{xu X2)) e(ti - t2) , (66) 

and 

T^xi)^X2))^ = D>{xi, X2) e(ti - t2) + D<{xi, X2) e(t2 - h) . (67) 

As discussed in Sect. 3, the time-ordered Green function is the sum of a 
disconnected and a connected piece, 

T^x,)Hx2))^ = {H^i))^ + ^A^^i' • (68) 
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a. 



Figure 6: Contributions to the self-energy 11 up to C(A^). 

In the retarded Green function the disconnected piece drops out. 

As usual the connected piece satisfies a Dyson-Schwinger equation, 

Dl{xi,X2) = Ga{xi-X2)+ j dyidy2Gtv{xi-yi)Tl{yi-y2)D\{y2,X2) . (69) 

The contributions to the self-energy are shown in Fig. ^ up to order A^. To 
leading order the self-energy is local. The contribution Fig. ^a is removed 
by a mass counter term and Fig. ^3 is the interaction with the background 
field, 

n(^)(2/i,2/2) = -i\ 5(2/1 - 2/2) . (70) 

The free propagator G\ may also be split into two contributions which 
describe propagation forward and backward in time, respectively, 

G'a(xi - X2) = G>(xi - X2) e(ti - t2) + G<(xi - X2) Q{t2 - h) . (71) 

From Eq. ( pOD and (|30|) one reads off, 

Gl{x,~X2) = / e-*(^-^^) + Gt,a(xi - X2) 

G<(xi -X2) = I e^'=-(-i--2) + Gt,a(xi - X2) . (72) 

The corresponding retarded Green function reads 

Df(xi-a;2) = -i{Gl{xi-X2-Gt{xi-X2))Q{ti-t2) 

= -2 I e^'^("^-"^) sin(a;k(ti - t2))Q{ti - . (73) 

Jk 
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Inserting Eqs. ( ]66| ) and ( [7I| ) in the Dyson-Schwinger equation ( |69| ) one 
obtains two coupled integral equations for D^{xi,X2) and D^{xi,X2). For 
the retarded Green function Dfi one finds, 

D^{xi - X2) = Gr{xi - X2) + ^\Jdy Gr{xi - y)^l{y)D^{y - X2) . (74) 

This equation is identical with the equation for the classical retarded Green 
function Hr obtained in ITDI. Instead of the classical field (j){y) now the 
expectation value 0a (2/) appears. 



a. 



b. 



Figure 7: Contributions to the classical self-energy 11 up to C(A^). 

In order to obtain the complete retarded Green function one has to per- 
form the integration over the soft modes, 



Dr{xi,X2) = J V(PVn e-f^"^ D^{xi - X2) 

uj<A 

This means, one has to evaluate the quantities 
J V(j)V'K e-^^^(f)l{xi) . 



(75) 



^'i{x2r, 



(76) 



This is completely analogous to the integration over initial conditions in 
For two fields one has to leading order in A (cf. Eq. 



V(I)Vtt e-^^*0A(a;i)0A(x2) oc GtA^i - X2) . 



(77) 
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For the products of 2n fields one obtains, 
J V<j)VTT e-'3^^0i(xi)...0i(x2„) 

a;<A 

= 2" ^{xi — X2) . . . ^{x2n-i — X2n) + permutations . (78) 

One easily verifies that the full retarded Green function again satisfies a 
Dyson-Schwinger equation, 

Dr{xi-X2) = Gr{xi-X2)+ J dyidy2 Gnixi-yi) Il{yi,y2) DR{y2, X2) . (79) 

The leading contributions to Tl{yi,y2) are shown in Fig. |^ The first con- 
tribution is subtracted by a mass counter term. The second contribution 
reads, 

miyi,y2) = \' Gniyi - 1/2) G|,^(l/i - 2/2) • (80) 

This yields 



(p, oj) = j dtd?x e'^'e-'^^ 11^^^ {t, x) 



1 1 



where r/j = ±1, i = 1, . . . , 3 and 



(81) 



P = ^3o ^3o ^3o 27r '5 p-qi-q2-q3 ■ 82 

[l-KyluJi [2'Ky2uJ2 [iTXyiuJ^ 

From Eq. (^) one obtains for the imaginary part of the self-energy (lc; > 0), 

r(2)(p,.;) = IX'T^ I d^{p)^^[uj,5{uj-u,-U2-oo,) 
/ J UJ1UJ2OJ3 

+ uj5{uj + uji — UJ2 — uj^)] . (83) 



As discussed in JlOl, this expression agrees with the imaginary part of the 



quantum self-energy to order for uj <^ T, where the Bose-Einstein distri- 
bution function becomes 

T 

n{iu) ~ - > 1 . (84) 



From Eq. ( p3D one easily obtains the plasmon damping rate, 

7 = — r(^) 0, m = 1 - — — + O — . 85 

' 2m ^ ' ^ 15367r m \ -k^ K \\K) J J ^ ^ 

Except for the corrections O (j^^ the result agrees with the quantum result 
to leading order in A The dependence on the infrared cutoff A controls 
the size of corrections to the infrared limit. 

To assess the size of the quantum corrections one has to compute the 
contribution of Fig. ^ to the damping rate (diagram c. does not contribute 
to the imaginary part, diagram e. containes less IR dominant propagators 
and therefore gives a smaller contribution). The correction to the classical 
self-energy n reads (cf. appendix) 

Sn{y^, y^) = y G^^iv, - y^) [G^+(yi - ^2) - G^Jz/i - ^2)] • (86) 

After some algebra one finds for the correction to the imaginary part of the 
on-shell self-energy (p = 0, w = m), 

51^2) (0, m) = — ttT / (i$(0) —5{m + U1-U2- uj^) 
2 J uj\ 

7687r 7r2 V 2A ^' ' J ^ ' 

At high temperatures ~ AT^, i.e., the quantum corrections are 0{-\fX\ii A). 
This implies that the classical limit gives a meaningful approximation to the 



quantum theory only to the leading order in A |10 



6 Summary 

We have shown that in scalar quantum field theories at high temperature 
the behaviour of low frequency modes is described by an effective classical 
theory. Corrections to this classical limit are controlled by the dependence 
of observables on an infrared cutoff A which is introduced to separate soft 
modes from hard modes. 

Integrating out the hard modes yields an effective theory which depends 
on the infrared cutoff A and on a classical background field (l)\{x), the expec- 
tation value of the field operator in a coherent state containing soft modes 
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only. The thermal Green functions of the low energy effective theory factor- 
ize in products of expectation values of the field operator up to corrections 
of order h(3K. This loss of quantum coherence is due to the interactions with 
the soft modes in the thermal bath. The behaviour of soft modes is described 
by classical statistical field theory. 

The separation of soft and hard modes has been carried out at some 
fixed initial time tj. Due to interactions with the thermal bath soft and 
hard modes will mix at later times. The implications for the behaviour of 
correlation functions at large times remain to be studied. 

Our analysis also illustrates some of the difficulties which have to be 
faced in calculations of real-time correlation functions in gauge theories at 
high temperature. The separation of soft and hard thermal modes is com- 
plicated by gauge invariance and in the cases of interest the gauge coupling 
is rather large, so that the magnetic mass scale g^T is not well separated 
from the temperature T. Hence, it is not clear to us how accurately existing 
classical numerical simulations for real-time correlation functions in gauge 
theories approximate the quantum theory. 

We would like to thank D. Boyanovsky, M. Liischer, A. Patkos and 
J. Polonyi for helpful discussions and comments. This work has been partially 
supported by the Hungarian NSF under contract OTKA-T22929. 



Appendix 

The contour integral in Eq. (^) is most easily evaluated in the two-component 
formalism which is obtained by taking the limits tj — oo, tf oo. Then 
only the contributions from the parts C+ and C_ of the Keldysh-type contour 
in Fig. 1 remain. 

The corresponding propagators read (cf. |[T6[| ) 

G++{Xi~X2) = e{h-t2)G>{Xi-X2) + Q{t2-h)G<{xi-X2) 

= (G__(xi-X2))* , (88) 
G+-{xi - X2) = G^{xi-X2) = {G -+{xi ~ X2))* , (89) 

where 

G>^<\x) = f ^^^(k,^) e'^" , (90) 
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= iG<{k,t)r . (91) 

The interactions of the two fields 0+ and 0_ which correspond to the two 
branches C+ and C_, respectively, are given by 

Si = -jdtjd^xj^({(t>+ + ct>!,Y-{ct>^ + ct>^f 



- I dt I d'xX i ^(0+ - + - <pl)<Pl 



+ l{<Pl'<P'-)<pA + ^{<PX-<Pt)) , (92) 

where 0a is the background field introduced in Sect. 2. 
For the contribution of diagramm |^ one now obtains 

r 

—(f)A{x) / dx' G{^{x,x'Y ^{x')\ 
4 J c 

^%A(a;) / rfV {G++{x - x'f - G+.{x - x'f) <Pl{x') 



4 
4 



<Pk{x) J dt'e{t - t') J dV [G>{x - x'f -G<{x- x')') (plix) . 



Inserting the expressions for the propagators G^ and and performing a 
partial integration with respect to t' one obtains the result Eq. (B3). 
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